On one integrable system with a cubic first integral 



A. V. Vershilov and A. V. Tsiganov 

St. Petersburg State University, St. Petersburg, Russia 
e-mail: alexander. vershilov&gmail. com , andrey. tsiganov@gmail. com 



PACS: 45.10.Na, 45.40.Cc 
MSC: 70H20; 70H06; 37K10 

f"^ ' Abstract 

Recently one integrable model with a cubic first integral of motion has been studied by Valent 
using some special coordinate system. We describe the bi-Hamiltonian structures and variables of 
separation for this system. 

1 Introduction. 
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The aim of this note is to consider one particular two-dimensional integrable model defined by natural 
Hamilton function 

2 

H = T+V= J2^(Q2)PiP j + V(q 1 ,q 2 ) (1.1) 
with the metric depending on one variable, and cubic additional integral of motion with the leading terms 



>: H 2 = pp\ + 2qT Pl + ... , pel, q>0. (1.2) 

According to [7J, relevant metrics are described by a finite number of parameters and lead to a large 
class of models mainly on the manifolds § 2 and H 2 . By suitable choices of the parameters entering 
the construction these systems are globally defined and contain as special cases the known systems of 
C*~) . Goryachev-Chaplygin, of Goryachev, and of Dullin and Matveev, see [7J . 

In 6] we introduce a natural Poisson bivector depending on two arbitrary functions, which allows us 
to describe similar family of integrable system with cubic additional integral of motion in spherical coor- 
dinates. Below we rewrite this bivector in the coordinate system used in [7J, calculate the corresponding 
variables of separation and prove that equations of motion are linearized on genus three non-hyperelliptic 
curve. 



2 Settings 

In this section we recall some necessary facts about natural bi-integrable systems on Riemannian manifolds 
admitting separation of variables in the Hamilton- Jacobi equation [3J [51 |S] 

Let Q be a n-dimensional Riemannian manifold. Its cotangent bundle T*Q is naturally endowed 
with canonical invertible Poisson bivector P, which has a standard form in fibered coordinates z = 
(q 1 , . . . ,q n ,p x , . . . ,p n ) on T*Q 

P=(_° I I), {/.«> = ^4f,*>=g^^, (2.D 

In order to calculate the variables of separation for the given integrable system with integrals of motion 
Hi , . . . , H n in involution 

{H i ,H j } = 0, i,j = l,...,n, 
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in the bi-Hamiltonian set-up we have to solve equations 



[P,P] = [P',P']=0, 
where [., .] means a Schouten bracket, and 

{i^l'-O, i,j = l,...,n, {/,.g}' 



(P'df,dg), 



(2.2) 



(2.3) 



with respect to the Poisson bivector P'. Then we have to calculate the so-called Nijenhuis operator (or 
hereditary, or recursion operator) 

P'P- 1 . 



N 



(2.4) 



If N has, at every point, the maximal number of different functionally independent eigenvalues ui, ■ ■ ■ , u n , 
they may be treated either as action variables (integrals of motion) or as variables of separation for this 
dynamical system [5J [5] . 

Separation of variables for natural integrable systems with higher order integrals of motion always 
involves generic canonical transformation of the whole phase space. The definition of the natural Hamil- 
tonians ([l.ljl is non-invariant with respect to such transformations of coordinates on the whole phase 
space (cotangent bundle). 

In the situation, when habitual objects (geodesies, metrics and potentials) lose their geometric sense, 
and the remaining invariant equations (|2.2I2.3[) have infinite number of solutions, the notion of the natural 
Poisson bivectors becames de-facto a very useful practical tool for the calculation of variables of separation 
[5J H] . It is an experimental fact supported by all the known constructions of the variables of separation 
on the sphere. We try to draw attention to this experimental fact in order to find a suitable geometric 
explanation of this phenomenon. 

Similar to the natural Hamilton function on T*Q, natural Poisson bivector P' is a sum of the 
geodesic Poisson bivector P' T and the potential Poisson bivector defined by a torsionless (1,1) tensor field 
A(gi, . . . , q n ) on Q associated with potential V [5J |B]: 



P' 



P't 



\ 







A* 



-A 



j * 



E 

fc=i 



<9A 



k'i 



dA 



dqj dqt 



Pk 



(2.5) 



The geodesic Poisson bivector P' T is defined bynxn geodesic matrix II on T*Q : 



P' 



fe=i 



x jfc(<z) 



an 



jk 



dpi 



Yik(q) 



dU lk 
d Pj 



V 



E 

k=l 



IT 



dn ki an k 



dqj dqi 



Zfc(p) 



(2.6) 



In fact, for the given matrix II functions x, y and z arc completely determined by the equations 

[P,P^] = [P^,P^]=0, (2.7) 

whereas A is obtained as a solution of the equation (|2.2[) . Discussion of this useful anzats (|2.5j) may be 
found in [51 [3] ■ 



2.1 Two-dimensional integrable systems with cubic integrals of motion 

Let Q be a 2-dimensional Riemannian manifold, and its cotangent bundle T*Q is naturally endowed with 
four canonical coordinates q\ t 2 and p\ 2. 
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One of the natural Poisson bivectors P' ([2. 5ft listed in [6] is defined by geodesic matrix depending 
on one variable q2'- 



I 



n : 



2 \dq 2 g(q2) J 



\ F 
and diagonal potential matrix 



F = (fl(«!)pa - i^fe)Pi) 2 , i = (2.8) 



Here <?(<?2) and ^.(92) are arbitrary functions and 




A = a exp ( i qi - I ^) dq 2 ) \ I , ael. (2.9) 



9(12) 

X22 



2ft(? 2 ) ' 

Other functions yik, ~x~ik and z k in the definition (12.61) equal zero. 

A characteristic polynomial of the corresponding recursion operator N = P'P^ 1 reads as 

det(N ~ 77I) = U 2 — (F + 2A U ) V + A-ii 2 



Here F (|2 .8[) is a complete square, it allows us to introduce the coordinates «i ;2 as roots of the linear in 
momenta pi.2 polynomial 

B(X) = (A — ui)(A - u 2 ) = A 2 - iVFA + An . (2.10) 

In this coordinates B = (77 + v 2 ) 2 (r] + v|) 2 . We prefer to use the linear in momenta polynomial B(X) 
instead of 8(77) because in this case it is easy to find a solution of the equations 

{B(X), A(m)} = 77— r f ^ - ^) , {A(X), A(m)} = , (2.11) 



/x — A \ A 

with respect to other linear in momenta polynomial 

idq2 ipi 
5(92) A 

Equations (|2.11|) entail that canonically conjugated momenta are equal to 



AW-/; 



p„ 1>2 - A(A = «i, 2 ). (2.12) 

Substituting variables 

x^agv^ 1 , z = \p k , fc = l,2, a eR, 
into the generic equation of the so-called (3,4) algebraic curve [T] 

$(z,x) = z 3 + (a 1 a;+a2)z 2 + (-ffi.T 2 +6 1 a;+fe2)^+a; 4 +-ff22; 3 +Cia; 2 +C2a;+c 3 =0, a k ,b k .c k GK, (2.13) 
and solving the resulting equations with respect to -ffi.2, one gets the following Hamilton function 

H X = T + V + J- ((a lW 2 2 - b lW2 + c 2 )h 2a lW2 + h) Pl + fefo"* ~ M» + C») (2 u) 

where geodesic Hamiltonian and potential are equal to 



T = 

alw 2 



■ — f (c 3 — 100 + a 2 w^ ~ b 2 w 2 )h 2 + (6 2 + 3w 2 ~ 2a 2 w 2 )h - 3w 2 + a 2 Vi 
,W9 \ ' ) 



[ V (ol 3, 2 . N, , r , n ,2 o ^ 5 2 (C3 - w| + 02^1 - 6 2 W 2 ) 2 

2(c 3 - w 2 + a 2 w 2 - o 2 W2)h + D2 + 3w 2 - 2a 2 W2 P1P2 5 P2 

V / aw? 



y _ alc~ lqi awi(c 3 - W2 + a2iv 2 - b 2 W2)e lqi c\ 

aW\W2 CIqW2 W2 
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Here 

( /%2) , \ f dq : 

Wi = exp - / — — - dq 2 , w 2 - ' 



5(92) / J g{q2) 

Second integral of motion H 2 is a cubic polynomial in momenta p\ i2 . 

If h and j ^ 0, in order to reduce Hamiltonian Hi (|2. 14[) to the natural form, we have to solve 
two integral equations 

a\w\ — b\w 2 + c 2 = , bi — 2a\w 2 = , (2-15) 

with respect to functions h(q 2 ), g{q 2 ) and parameters a±,bi,c 2 . If we want to obtain diagonal metric, 
then we have to add one more equation to this system. There is an additional freedom related with a 
possible canonical transformation p 2 — > p 2 + / (q 2 ) , which change linear in momenta terms in (|2.14[) . 
If Q = S 2 is a two-dimensional unit sphere with spherical coordinates 

9 = (91,92) = 0,#) and P=(pi,P2) = (p<t>,Pe), (2.16) 

then we could get a whole family of natural integrable systems on the sphere using different functions 
h(9), g(8) labeled by particular values of parameters in (|2.15[) [HIIH]- 

3 Natural Poisson bivectors in (^-variables 

In [7j, solving the equation 

{H u H 2 } = 

in framework of the Laplace method, author prefers to use special canonical coordinates 

9i=0, 92 = C, Pi=P<p, P2=Pc, (3.17) 



Here Hi is a natural Hamilton function (II. ip with diagonal metric g(£) and H 2 is fixed by fL 

According to [7], special case q = in (II. 2\i is rather difficult to obtain as the limit of the general 
case 9 7^ 0. So, we first integrate the special case at q = and only then we consider the generic case. In 
contrast with [7], we restrict ourselves only by local analysis. 

3.1 Case q = 

Let us take integrals of motion from 7 , Theorem 1 formulae (14-15): 



2 V " 

H^ 0) = pl-2X (VF sin # f + (VF)' cos (p P4 ) - 2 m , 



(3.18) 



with 

F = -2p + 3c C + C 3 , G = 9c 2 + 24p C - 18c C 2 - 3C 4 • (3.19) 

Here A, /z, po and Co are arbitrary parameters. 

Substituting new variables (|3.17p into definitions (|2.8I2.9[) one gets a bivector P' depending on two 
arbitrary functions h(C) and 5(C)- Using this bivector we can easily solve the equation (|2.3[) and prove 
the following proposition. 

Proposition 1 Integrals of motion h[° 2 \3.18)) are in bi-involution 

} = {<Uf}' = (3.20) 
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with respect to a pair of compatible Poisson brackets associated with the canonical Poisson bivector P 
\2. 1\) and natural Poisson bivector P' \2. 5\) defined by 



n 



d_ 



EL 

~F 



V o 



and 



F 

4Ae-^ 



1 
1 



iF' 



F 

X22 = -JJ- 



(3.21) 



(3.22) 



As a consequence, substituting = 2 and h(() = (InF)' into variables Wi j2 (|2.10p and p vi 2 (|2.12[) one 
gets variables of separation for this integrable model. 

Below, in order to extend the known palette of natural Poisson bivectors listed in [B], we consider 
another solution of the equation (|2.2j) depending on both variables 4> and £. 



Proposition 2 Integrals of motion h[°2 13.18\) are in the bi-involution &3.20\) with respect to canonical 

Poisson bracket and bracket {.,}' associated with the natural Poisson bivector P' 12. 5\) defined by 2 x 2 
geodesic matrix 

iF' \ 

F = 



n 



F 



iF' 



2^F 



P4> 



FPC 



diagonal potential matrix 



and function 



A = Xc^VF ( J J 



(3.23) 



(3.24) 



yn 



Other functions y^, and Zfc in the definition 



equal zero. 



As for the Lagrange top, this geodesic matrix (|3.23j) is factorized LT = e 2l< ^n(C), see discussion in [6j. 

We are now able to analyze the corresponding recursion operator N = P' P~ x . For instance, we can 
define the variables of separation u\ 2 as roots of the linear in momenta polynomial 



B(v) 



(77 - ui)(r) - u 2 ) — if — h/F r] + Ai 



(3.25) 



ry 2 - ie* 



n / ! l(CO+ o C2 ? ^P0 + V / -2po + 3coC + C 3 PcU + Av / -2po + 3coC + C 3 e i0 , 
2^-2po + 3c C + C 3 / 



so that the characteristic polynomial of N reads as det(N — 77I) = (77 — u\) 2 (jj — u 2 ,) 2 . The conjugated 
momenta are equal to 

_ i^_iAC 

The inverse transformation looks like 

\u\U2(uip Ul - "2P« 2 ) _ K u iPu! - u\p U2 ) 



C 



X[U\ — U 2 ) Ui — U2 

(3.26) 



i, /A 2 F\ iF' iX( Ul +u 2 ) 

2 \U\U2) 2F wiu 2 
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In these variables of separation our initial bivector P' (13.21|3.22j) looks like 

/ 

P' =3 

V 

Now we can calculate the integrals of motion 

(0) _ \pu\u 2 (uip Ul ~u 2 p U2 ) (u 2 + uiu 2 + uppoX 2 3ic (uip U2 - p Ul u 2 )X niit 2 









2 2 


\ 











2 2 


2 2 



















o 1 



H, 



(0) 



A(ui — u 2 ) u iU 2 2(«i — 1*2) 

K^iPui +u 2 p U2 uip Ul + ulpl 2 )(uip Ul -u 2 p U2 )uiu 2 
2(ui — u 2 )X 

2\p(p Ul u\ - p U2 u\) 2(u 1 +u 2 )p X 3 3ic (p Ul -p„ 2 )A 2 



ui — u 2 
and the control matrix 

U\U 2 



l l u 2 



Ml — U 2 



X 




from the definition P'dH\ 0) = P ^ Fy dff 

3=1 



(0) 



,2. 



Suitable normalized left eigenvectors of control matrix F form the Stackel matrix S, so the notion of F 
allows us to compute the corresponding separated relations 



i=i 



1,2. 



Here functions and Ui depend only on one pair (ui,p Ui ) of canonical variables of separation [3J. 

In our case the integrals of motion and the variables of separation are related via the following 
separated relations 



z) = z 3 + (3A 2 c - 2pu 2 )z + Xu 4 - H { 2 °'u 3 + 2XH{>u 2 - 2X 3 p a = 0, 



(3.27) 



at u = ui. 2 and z = iu\ 2 Pi, 2 - Equation z) = defines genus three non-hyperelliptic curve with the 
following base of the holomorphic differentials 



du 



2/j.u 2 — 3A 2 co — 3z 2 ' 



n 2 



udu 



2jj,u 2 — 3A 2 cq — 2>z 2 ' 



n 3 = 



zdu 



2p,u 2 - 3A 2 c - 3z 2 ■ 



According to [T] it is the so-called (3,4) trigonal curve. It is easy to prove that the equations of motion 
in variables of separation have the following form 



Ul 



U-2 



2pu\ - 3A 2 c - 3zj 2pu\ - 3A 2 c - 32 



— , z k = m k p k 



U\U\ 



u 2 u 2 



2p,u\ - 3A 2 c - 3z 2 2/mf. - 3A 2 c - 3z| 
The aforementioned quadratures in the integral form 



fi 2 



08, 



(3.28) 



represent the Abel-Jacobi map associated to the genus three non-hyperelliptic curve defined by the 
equation $(q,p) = (|3.27|) . In particular it means that instead of z in O12 (j3.28[) we have to substitute 
the function on u obtained from the separated relation p.27[) . 



G 



3.2 Case q 

Let us take integrals of motion from 7 , Theorem 5 formulae (39-40): 



Hi = -(fpI + — pSV— W + A, 
#2 = PP^ + 2qH lP(t , - \/Fsin0p c - (Vf)' cos 



(3.29) 



with 

F = C o + Cl C + c 2 C 2 + C3C 3 , G = F' 2 — 2FF" , c 3 = ^ ( 3 - 3 °) 
where (3q, Co, ci, c 2 and p, g are arbitrary parameters. 

Proposition 3 .A/ter substitution of the new function F $3. 30\) into the previous definitions \3.21\3. 22\) 
and i3.23\3.2~4\ ) one gets two compatible Poisson bivectors P' and P' , so that the integrals of motion .Hi, 2 
H3.29\) are in the involution with respect to the corresponding Poisson brackets. 

Other calculations are standard. Namely, the variables of separation are given by 

— p lp + VFp c )n-?L_, (3.31) 

iP0 C_ 

u k 2qu% 

whereas the inverse transformation reads as 

_ 2quiu 2 {uip Ul - u 2 p U2 ) „ _ K u iPu! ~ ulp U2 ) 



U\ — U2 U\ — U2 

1 , / Aq 2 u\ul\ IF' u\ + u 2 

2 \ b J 2b 2qu\U2 



(3.32) 



Then we have to calculate the control matrix F and the corresponding Stackel matrix S in order to get 
the desired separated relations, which are obtained by substituting u = ui t 2 an d z — pu\ 2 p Ul 2 m the 
following equation 



^ O 3 2 (8q 2 H lU 2 - C]> 2 4 2 3 /3 P 2 M 2 J3 2 C 

<±> = 2z — C2 z z — p u — 2\p bi-2,u 7t = . (3.33) 

2q ' q <iq z 



This equation z) = defines the genus three non-hyperelliptic curve ( so-called (3,4)-curve, see pQ) 
with the holomorphic differentials 

Ul = ~2c^/dz ' " 2= -2c 3 9$/5z' " 3 = -2e 3 d*/dz ' ( } 

The desired quadratures read as 

M^l UU2 _„ _ 2 

8p 2 H lU 2 - 12c 3 z 2 + 4C2C3Z! - cic§ 8p 2 ffi U 2 _ + 4c 2 c 3 Z2 - cic 2 ' z >°-P u kPk, 

and 

ZlUl 22^2 P 



8p 2 H\u\ — 12c 3 z 2 + 4c2C 3 zi — cic§ %p 2 H\u\ — 12c 3 z| + 4c 2 c 3 z 2 — cic 2 8 ' 
The above quadratures in the integral form 

r> 2 + / n 2 = A, / n 3 + n 3 = -t t + p 2 , (3.35) 
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represent the Abel- Jacobi map associated to the genus three non-hyperelliptic curve defined by (|3.33[) . If 
we change u —> iu, that corresponds to transformation P' —> —P', one gets the equation <E>(u, z) (|3.33[) 
with real coefficients, but the coefficient before time variable in f|3 . 28[) becomes imaginary number. 

In order to give an explicit theta-functions solution, one could apply the standard Weierstrass ma- 
chinery describing the inversion of the hyperelliptic quadratures [10 . However, for the genus three 
trigonal curve such solution of the Jacobi inversion problem is only created. Some of the results referring 
to Weierstrass theory for the general trigonal curve of genus three may be found in pQ. 

On the other hand, for a non-hyperelliptic Riemann surface of genus three the moduli spaces of rank 
two bundles compactifies into a singular projective variety, which is closely related to the Kummer variety 
and which is a quartic hypersurface of P 7 . We can try to apply the generic theory of vector bundles on 
algebraic varieties to investigation of this integrable system and vise versa [8] . 



3.3 Change of the time 

If the dimension of the configurational space n is less then genus g of the separated curve, then the Abel 
map is either lack of uniqueness or degenerate. In this case we are free to choose n integrals of motion 
from the coefficients of the separated curve, that equivalent to choose the corresponding time variables. 
The Kepler change of the time, the Liouville transformations, the Maupertuis- Jacobi transformations 
and the coupling constant metamorphosis are examples of such duality. One of the well-known examples 
is a duality between harmonic oscillator and the Kepler system, which was expanded on the quantum 
case by Schrodinger. 

All the known such transformations are associated only with the hyperelliptic curves, see [3J[3]. Here 
we want to show one example related with non-hyperelliptic curve. Let us consider the following equation 

* 1 3 2 (8q 2 au 2 ~ d)p 2 4 25 3 2 - 2 p 2 C , . 

<P = 2z — C2 z z — p u — zip H2U + 2p H\u — = , (3.36) 

which is related with the initial equation (|3.33|) by permutation of coefficients a and H\ only. If we 
substitute variables (|3.31j) into ()3.36|) and solve the resulting pair of separated equations, then one gets 
the Hamilton function 

F = c + ciC + c 2 C 2 + c 3 C 3 , G = F' 2 -2FF'. (3.37) 



Second integral of motion in this case is equal to 



3 F' /— 2ap 

H 2 = PPd - — 7= P<p cos cp - VFp c sin 4 ' 

ZV r 



C3 

Integrals fli j may be obtained from (|3 . 1 8|) at special values of parameters using canonical trans- 
formation £ — > £ — C2/3C3. 

Dual Stackel systems with Hamiltonians H\ and Hi are related to each other by a canonical trans- 
formation of the extended phase space [21 H] ■ For the both systems trigonal curve defined by (I3.33|) or 
(|3.36p has the same holomorphic differentials fij. (|3.34[) . but instead of quadratures (|3.35|) in the second 
case we have 

/•Ml rU 2 1 /"ill f«2 

/ fix + / n x = - — t + fa , / n 2 + o 2 = ft , (3.38) 

Similar to hyperelliptic case, this duality (change of the time) is related with the different choice of the 
holomorphic differentials in the inversion Jacobi problem [31 13] . 

Using the same permutation of coefficients in the generic equation (|2.13l) of (3,4)-curve we can obtain 
the dual Hamiltonian to more complicated Hamilton function f|2 . 14|) . Other possible generalization is 
related with an another pair of differentials in the Abel map: 

1. case q — is associated with ft 1 and p.38[) ; 
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2. case q ^ is associated with O2 and p.35j) : 



3. one more integrable case is associated with fix and O3. 

Some difficulty here is related with the coordinates choice, because we can not directly recognize inter- 
esting physical models in ^-variables. So, different possible generalizations of the Goryachev, Chaplygin, 
Dullin-Matveev systems and other systems from [5J|H] associated with the generic (3,4)-curve in term of 
physical variables on the two-dimensional sphere will be discussed in a forthcoming publication. 

According to [6] , there are natural generalizations of matrices II and A (|2.8|2.9j) on the 3-dimensional 
case. It can allow as to obtain a new three-dimensional natural integrable system with higher order 
integrals of motion constructively. 
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